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ABSTRACT 



Aims. We study the dynamics and the evolution of a C2.3 two-ribbon flare, developed on 2002 August 11, during the impulsive 
phase as well as during the long gradual phase. To this end we obtained multiwavelength observations using the CDS 
spectrometer aboard SOHO, facilities at the National Solar Observatory/Sacramento Peak, and the TRACE and RHESSI 
spacecrafts. 

Methods. CDS spectroheliograms in the Fe xix, Fe xvi, O v and He I lines allows us to determine the velocity field at different 
heights/temperatures during the flare and to compare them with the chromospheric velocity fields deduced from Ha image 
differences. TRACE images in the 17.1 nm band greatly help in determining the morphology and the evolution of the flaring 
structures. 

Results. During the impulsive phase a strong blue-shifted Fe xix component (—200 km s~ ) is observed at the footpoints of the 
flaring loop system, together with a red-shifted emission of O v and He I lines (20 km s _1 ). In one footpoint simultaneous Ha 
data are also available and we find, at the same time and location, downflows with an inferred velocity between 4 and 10 km s -1 . 
We also verify that the "instantaneous" momenta of the oppositely directed flows detected in Fe xix and Ha are equal within 
one order of magnitude. These signatures are in general agreement with the scenario of explosive chromospheric evaporation. 
Combining RHESSI and CDS data after the coronal upflows have ceased, we prove that, independently from the filling factor, 
an essential contribution to the density of the post-flare loop system is supplied from evaporated chromospheric material. 
Finally, we consider the cooling of this loop system, that becomes successively visible in progressively colder signatures during 
the gradual phase. We show that the observed cooling behaviour can be obtained assuming a coronal filling factor of « 0.2 to 0.5. 
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1. Introduction 

Most solar flare models are based on the idea that flares 
arise from the sudden release of free magnetic energy 
stored in the corona i n non-potential configurations (e.g., 



Klimchuk et alJ ll988). The bulk of the energy is released 



through reconnection of the magnetic field lines in the 
coron a ( PetschekHl964t iKopp fe Pneuman 119761: IShibatal 
19961) and then transported down to the chromosphere 
along the magnetic field lines by accelerated particles or 
by a thermal conduction front. Hydrodynami c simulations 
indica te that both non-ther mal 
1985aflbl ; [Mariska et alj Il989h 



electrons (jFisher et al 
and a thermal conduction 
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front (|Gan et al.lll991l ) should lead to the evaporation o f 
chromospheric material. Models by iFisher et al. (|l985a[) 
result in upflows (gentle evaporation) or downflows (ex- 
plosive evaporation), on both transition-region (TR) and 
chromosphere, depending on the energy flux of the imping- 
ing electrons, with downflows predicted for large energy 
fluxes. In the case of explosive evaporation the overpres- 
sure of evaporated plasma drives upward motions in the 
corona and downward motions in the chromosphere and 
it is expected that the momenta of the oppositely mov- 
ing plasma should be balanced. The chromospheric evap- 
oration injects mass into the corona to fill coronal loops 
that appear as the soft X-ray (SXR) flare and that even- 
tually cool down radiatively and/or conductively form - 
ing cold Ha post-flare loops (see, e.g.. lMoore et al. I ll980h . 
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Although this general picture seems well established, the 
details of some of the key processes assumed to take place 
remain to be verified. The small temporal and spatial 
scales involved, together with the need for comprehensive 
temporal and height coverage, conjure up to make "com- 
plete" flare observations a very demanding task. 

Blue-shifts of lines formed at flare temperatures dur- 
ing the impulsive phase (corresponding to upflow ve- 
locities in excess of several hundreds kilometre per sec- 
ond) have been observed since the ear ly 1980s (e.g., 
Doschek et alll980tlAn"tonucci et al.lll982f) and confirmed 
with t he Bragg Cry s tal S pectrometer aboard Yohkoh 



(BCS, ICulhane et al. 199ll ). However, the majority of 



these observations revealed only a blue asymmetry, in- 
dicating the presence of a large, static co mponent even 
in the early phases of the flare (e.g., iBentlev et al 



1994 iGan fc Watanabel[l997h . rather than wholly blue- 
shifted profiles as predicted by hydro dynamic models (e.g. 



Emslie et alJll992t IGan et al.lll995f ). To which extent the 



lack of spatial resolution might have influenced these ob- 
servations is not easily quantifiable (neither BCS nor ear- 
lier spectrometers provided any effective spatial resolution 
on the solar surface). Since reconnection is expected to 
occur at very small spatial scales, it might well be that 
the associated chromospheric evaporation takes place in 
rather confined portions of the solar surface. This is indeed 
the case in eruptive flares, for which, at any given time, 
the evaporation in the reconnection scenario occurs in a 
very thin shell ( < 1") located on the outer part of the flar - 
ing loop system (jForbes fc Actonlll996l : lFalchi et al.lll997l ). 
and the heating of i ndividual field lines is supposed to 
last only 100 s or so (|Forbesl 1200.4 The smearing intro- 
duced by observing over the whole disk would then signif- 
icantly alter the relevant signals. Hydrodynamic simula- 
tions mode l ling solar flare emission with m ultiple loops 
(jHori et all [l998t IWarren fc Doschekl l2005h rave shown 
that it is possible to reproduce the line profiles dominated 
by the stationary component (as observed by BCS) with 
a simple model based on the successive independent heat- 
ing of small scale threads. Only the simulated line profiles 
for individual threads do show strongly blue shifted line 
profiles. 

Spatially and temporally resolved observations of flare 
lines have finally become available with th e advent of the 
Coro nal Diagnostic Spectrometer (CDS, IHarrison et al 



19951) aboard SOHO, which provides stigmatic slit im- 
ages in the 31 to 38 nm (NIS-1) and 51 to 63 nm (NIS- 
2) range with a spatial resolution of few seconds of arc 
and a cadence depending upon the adopted observing pro- 
gram. Various authors have since reported the presence of 
blueshifted emission in hot coronal line s during flares, sig- 
nalling upflows of up to ~ 200 km s" zaykowska et al 
1999tlDel Zanna et al l 120021: iBrosiu sl 120031; iTe riaca et al 



events analysed up to now remains small. Moreover, the 
same authors report about intricate patterns of blue- and 
red-shifted emission in TR lines, depending both on the 
location within the flaring region and the phase of the 
flare, that still need to be consistently explained within 
the theory of chromospheric evaporation. Finally, spatially 
and temporally resolved observations of flare lines need 
to be combined with simultaneous chromospheric obser- 
vations, in order to verify the spatial and temporal co- 
incidence of chromospheric downflows with hot plasma 
upflows. In fact, while many authors have reported on 
the simultaneous appearance of redshifts in chromospheric 
lines and blue- shifts in flare lines dur i ng the early phase 
of flares (e.g.. ICanfield et alj Il990bl : iFalchi et alj Il992 : 
Wiilser et al. 1994) . to date the only direct proof that such 
shifts ari se from spatially co incident areas has been pro- 



vided bv lTeriaca et al.l (|2003l ). However, due to the obser- 
vational setup, their coronal and chromospheric data were 
not simultaneous, leading to the need of further observa- 
tions. 

In this context, we report here on comprehensive, 
multi-instrument observations of a small eruptive flare, 
covering its temporal evolution from the impulsive phase 
and related evaporation throughout the gradual phase 
and appearance of cold postflare loops. Spatially resolved 
CDS spectra allow us to study the dynamics of the up- 
per atmosphere at various times during the flare, while 
simultaneous chromospheric observations provide a test 
for the chromospheric evaporation scenario. RHESSI data 
acquired during the gradual phase, combined with CDS 
measurements, provide a final element in support of the 
chromospheric evaporation. Other SOHO instruments and 
TRACE, provide various context information and help us 
in defining and constraining the flare evolution. 



2003$ iBrosius fc Phillips! |2004) , but the total number of 



2. Observations and data analysis 

The observations were acquired during a coordinated 
campaign between ground based and SOHO instruments 
aimed at studying flare events by sampling the solar 
atmosphere from the chromosphere to the corona. The 
main instruments involved were the CDS/NIS spectrom- 
eter aboard SOHO and facilities at the NSO/Sacramento 
Peak. The analysed flare developed at about 14:40 UTC 
in region NOAA 10061 on 2002 August 11 (N10W20). 
Although small (GOES class C2.3), the flare was clearly 
eruptive in character with an impulsive rise and a long 
decay phase, and developed two ribbons visible at vari- 
ous wavelengths on opposite sides of the magnetic neutral 
line. A second flare, not studied here, developed in the 
same region around 16:25 UTC. 
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Fig. 1. TRACE image (logarithmic scale) in the 17.1 nm 
band showing AR 10061 at flare peak. The dark grey (red 
in the electronic version) box indicates the UBF FOV 
while the area rastered by CDS/NIS is represented by the 
light grey (blue) box. The white vertical line marks the 
position of the DST spectrograph slit. Here and in the fol- 
lowing images, the white dashed line shows the apparent 
magnetic inversion line obtained from MDI data in the 
region around the filament. 

2.1. Data 

Ground-based data - Monochromatic images at several 
wavelengths have been acquired, with a temporal cadence 
of a few seconds, by the tunable Universal Birefringent 
Filter (UBF: Ha line centre, —60 pm, +60 pm and 
— 150 pm off centre; He I D3; Na I D2 and in the con- 
tinuum) and by the Zeiss filter (Ha +150 pm) at the 
Dunn Solar Telescope of the National Solar Observatory 
(DST/NSO). The field of view (FOV) was of about 
150" x 150" with a spatial scale of 0.5" x 0.5". Spectra have 
been acquired with the Horizontal Spectrograph (HSG) in 
three chromospheric lines (Ca II K, H7 and He 1 D3) with 
a temporal cadence of 4 s and the fixed slit positioned as 
indicated in Fig. [1] 

Space-born data - Spectroheliograms of a 120" x 110" 
area (mostly overlapping the UBF FOV, see Fig. []} were 
obtained with the Norm al Incidence Spectro meter (NIS) 
of the CDS experiment ([Harrison et al. 19951 ) starting at 
14:26:27 UTC. Each raster was obtained in « 210 s by 
stepping the 4" wide slit eastward in 19 steps of 6". On 
board binning over two rows yielded a pixel size of 3.4" 
along the slit. Spectra were acquired in the He 1 58.43 nm 



(2 x 10 4 K), O v 62.97 nm (2.5 x 10 5 K), Fe xvi 36.08 nm 
(2 x 10 6 K) and Fe xix 59.22 nm (8 x 10 6 K) lines. Every 
ten rasters an additional step of ss 5" westward was per- 
formed to compensate for the solar rotation. During the 
entire period of CDS observations AR 10061 was con- 
stantly monitored by the Transition R egion and Coronal 
Explorer (TRACE. Itlandv et al.lll999h that provided im- 
ages in the 17.1 nm band (Fe IX - X, « 10 6 K) with a 
spatial scale of 1" pixel -1 (on board 2x2 binning) and 
a cadence of ss 7.5 s. The Reuven R amaty High-En ergy 
Solar Spectroscopic Imager (RHESSI, iLin et al.ll2002l ) ex- 
ited its night a few minutes after the peak phase of the 
flare, and provided information on hot (« 10 7 K) ther- 
mal flare plasma as well as on non-thermal components. 
From RHESSI data in the energy band 3 to 10 keV, we de- 
rived images with spatial resolution up to 2" and spectra 
of the spatially integrated flux with 1 keV spectral resolu- 
tion. Unfortunately, due to particle precipitation events, 
only three different times could be analy sed during the 



gradu al phase of the flare. Finally, MDI (jScherrer et al 



I1995I) full-disk magnetograms taken at 14:27:34 UTC and 
16:03:34 UTC were used to determine the position of the 
photosphcric apparent magnetic inversion line. 

2.2. Data co-alignment 

Data from the different instruments were aligned using 
MDI continuum data as a reference. UBF continuum im- 
ages (and, hence, all DST data) were easily aligned first. 
Afterwards, CDS images in the He 1 58.4 nm line were 
aligned with the UBF images acquired in Ha line cen- 
tre. Moreover, the flare kernels at the loop footpoints in 
the TRACE 17.1 nm images appear strikingly similar to 
those observed in Ha line centre, allowing a very good 
relative alignment. The brightening of the flare kernels is 
also very clear in the CDS O v images, providing further 
constraints. We estimate the alignment between TRACE, 
MDI and ground data to be precise within 1" and within 
2" with CDS. RHESSI images are more difficult to align 
and the precision is estimated around 5". Each spectrum 
or image has been tagged with its precise acquisition time 
so that light and velocity curves obtained using data from 
different instruments can be properly computed. 

2.3. Measurement of line-of-sight motions 

Ground-based data - Strong flows in the flaring chromo- 
sphere are often measured through spectral line asym- 



metr i es, defined general ly with the bisector (jZarro et al 



119881 : iFalchi et alJ 119921 ). In the portion of the flare rib- 
bon intersected by the HSG slit, the Ca 11 K line shows 
indeed an emission core with a red asymmetry. The H7 
line, instead, remains in absorption (with a central ra- 
diance higher than in the reference quiet area), so that 
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several blends in its wings prevent a reliable measure of 
the otherwise visible asymmetry. The He I D3 line is only 
barely visible. The velocity component along the line of 
sight (LOS) for the HSG spectra was hence computed us- 
ing the bisector of the emission core of the Ca n K line, in 
all 16 pixels in the ribbon area and for 41 different times. 
Unfortunately, the (fixed) spectrograph slit was positioned 
slightly outside the FOV covered by the CDS rasters (see 
Fig. [T|), so that no direct comparison of velocities deduced 
from HSG spectra can be performed. However, in the area 
of the flaring ribbons within the CDS FOV, we could use 
the Ha images acquired at several wavelengths to esti- 
mate the chromospheric velocities. To this end, difference 
images at ±60 pm and ±150 pm were computed. In all 
flaring kernels the differences between the red and blue 
wings of Ha were positive and very close in value, while 
the radiances in both the centre and in the wings remained 
higher than in the reference quiet area. This implies an 
asymmetric Ha profile with a stronger emission in the 
red part of the line and hence a downward velocity for 
all flaring kernels. To estimate the value of this velocity, 
we calibrated the radiance differences obtained from the 
Ha images in the kernel where HSG spectra have been 
acquired with the velocities obtained from the Ca II K 
spectra. This was done for all times and pixels for which 
a reliable estimate of the Ca n K velocity was available 
(about 500 points). We find that the average Ha radiance 
difference observed in the flare kernels, both at AA±60 pm 
and AA± 150 pm, corresponds to a downward velocity be- 
tween 4 and 10 km s _1 . Similar velocity values have been 
found from Ha line pro files in the case of sm all (GOES 
class B - C) flares (e.g., Schmieder et al.lll998 ). 

CDS spectra - Spectra acquired after recovery of the 
SOHO spacecraft are characterised by broad and asym- 
metric line profiles that were fitted with an opportune 
template provided within the CDS softwar^]. The tem- 
plate consists of a Gaussian component (with all param- 
eters free to vary) and a wing (non-Gaussian) component 
whose parameters are established for both NIS-1 and NIS- 
2 spectra. 

The analysis of line positions reveals that no trends 
(either parallel or perpendicular to the slit direction) are 
present within a single raster. However, the analysis of the 
raster-averaged NIS-2 line positions versus time reveal a 
small linear trend that is equal in He I and O v. The 
shift amounts to w 20 km s _1 over two hours. This trend 
has been accounted for by providing a corrected (shifted) 
wavelength vector for each NIS-2 line in each raster. No 
such a trend is visible in the Fe xvi NIS-1 line. 

Flows in the He I 58.43 nm, O v 62.97 nm, and 
Fe xvi 36.08 nm lines have been measured fitting a single 
component to the line profiles. These lines have an obvi- 
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Fig. 2. Soft X-ray full disk flux at 1 UA in the 
1 to 8 A band around the time of the flare (GOES 8 data). 
The time derivative (dashed line) represents a proxy for 
the hard X-ray curve (Neupert effect) and shows its max- 
imum around 14:41:20 UTC. 

ous pre- flare component, so that a reference wavelength 
has been obtained averaging the central positions derived 
from the fitting for the whole dataset. Velocities have then 
been computed simply using the Doppler-shift of the flare 
profiles with respect to the reference wavelength. 

The measurement of flows in the hot (log T/K=6.9) 
Fe Xix 59.22 nm line requires a more careful analysis. 
The line appears only during flares, so great caution is 
ne eded in choosing the reference p rofile (see discussions 



Teriaca et aLll2003tlBrosiusll2003l ). The line profiles ob 



1 CDS Software Note 53 by W.T. Thompso n 
http : // solar . bnsc . rl . ac . uk/ software/notes . shtml 
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tained averaging all CDS pixels with a radiance greater 
than 0.16 W m~ 2 sr _1 are identical in width and position 
in three successive scans during the late phase of the flare 
(from 14:53:10 to 15:00:10 UTC, times at raster centres), 
so their average was used as our reference profile. A multi- 
component fitting was hence applied to the profiles in var- 
ious pixels and times, constraining the rest component 
by imposing width and position of the reference profile, 
while the ba ckground (characterised by the Fe XII 59.26 
nm line, see iDel Zanna fe Mason I [2005) was constrained 
through the average non-flare spectrum. Uncertainties on 
best-fit parameters are evaluated accounting for data noise 
(photon statistics, pulse height distribution and read- 
out noise) and the average 3-cr uncertainty amounts to 
15 to 18 km s _1 over the active region for single compo- 
nent fitting and to 25 to 30 km s _1 for double component 
fitting. 

3. Flare evolution 

Since RHESSI was in the Earth shadow during the initial 
phase, no hard X-ray observations were available to define 
the impulsive phase of the flare. Moreover, the low tern- 
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Fig. 3. Evolution of the flare: TRACE 
images (logarithmic scale) in the 
17.1 nm band are shown in the first 
two rows and Ha images in the last two 
rows. The white box in TRACE images 
indicate the Ha FOV. The images at 
14:34 UTC show the region just before 
the flare. The areas (labelled A and B) 
over which data were integrated to ob- 
tain the light and velocity curves shown 
in Fig. |4] are indicated in black on the 
TRACE image taken at 14:42:00 UTC. 
The arrow in the last Ha image indi- 
cates the cool post-flare loops. 



poral resolution attained with the CDS rastering mode 
prevents the use of EUV radiance peaks as a p roxy for 
hard X -ray (HXR) bursts, as was done, e.g., bv iBrosiua 



(2003). Hence we used the time-derivative of the GOES 
SXR flux as a proxy for the HXR emission (Neupert effect; 
Neupertl[l968h . Both the 0.1 to 0.8 nm (1.6 to 12.4 keV) 
flux and its derivative are displayed around the time of 
the flare development in Figure [2] The main episode of 
energy release indicated by the Neupert peak is around 
14:41:20 UTC, but significant signal is present up to the 
GOES maximum, i.e., up to around 14:48 UTC. We as- 
sume the gradual phase starts after this time. 

Figure [3] outlines the main phases of the flare evolution 
as seen in the TRACE 17.1 nm channel and Ha centre. 
A filament located on the magnetic neutral line is clearly 



visible in both spectral ranges before the flare. The ac- 
tivation of a portion of this filament is observed around 
14:13 UTC in the Ha images and lasts up to 14:37 UTC, 
when it becomes clearly visible also in the TRACE im- 
ages. The eruption was detected by CDS in the O v line, 
where upward speeds around 30 km s _1 were detected at 
14:36:33 UTC. 

Few minutes later the footpoints of a large loop sys- 
tem (spanning over 6 x 10 4 km on the solar surface) start 
to brighten, reaching their maximum around 14:42 UTC. 
The flare ribbons are clearly visible and nearly identical 
in both TRACE and Ha images. Although the east-most 
ribbon is outside of the FOV of the UBF Ha images, we 
verified its existence in th e images acquired by th e HASTA 
telescope (see Fig. 1 in Maltagliati et al. 2006). During 
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Fig. 4. Light and velocity curves at the 
footpoints (areas A and B in Fig. [3} of 
the loop system. For clarity reasons un- 
certainties are only shown for the Fe xix 
speeds. O V, Fe xvi, and TRACE radi- 
ances are normalised to pre-flare values 
while Fe xix radiances are normalised 
to the average value reached after the 
flare (around 15:32 UTC). The aver- 
age uncertainties in the normalised ra- 
diance are of about 0.15. On footpoint 
A, the average uncertainties in speed 
are 12 km s" 1 (O v) and 20 km s" 1 
(Fe xvi). On footpoint B, the aver- 
age uncertainties in speed are 9 and 
12 km s" 1 for O v and Fe xvi. The 
shaded area represents the time dur- 
ing which downflows between 4 and 
10 km s" 1 were observed in Ha. (We re- 
mind that only for footpoint B we have 
chromospheric observations). On both 
panels, the scale of downward (positive) 
velocities has been expanded for sake of 
clarity. 



the impulsive phase of the flare, the ribbons are associ- 
ated with large coronal upflows and chromospheric and 
TR downflows (see next section). The brightenings last 
several minutes and fade off around 14:50 UTC in the 
17.1 nm band while their decay is much longer in Ha. 

Immediately after the GOES maximum, the first 
RHESSI spectrum (14:50 UTC) between 4 and 10 keV 
shows a non-thermal component indicating that electron 
acceleration is still present. During the gradual phase of 
the flare, two cooling loop systems connecting the ribbons 
appear in the TRACE images, one around 15:40 UTC, 
and the other around 16:03 UTC. Fig. [3] shows them at 
the time of maximum radiance. Finally, these loops appear 
also as dark features in the Ha images at 16:21 UTC, re- 
vealing that the plasma has cooled down to temperatures 
around 10 4 K. 

In the following, we will concentrate on the footpoints 
of the second loop system, contained in the CDS FOV, 
where most of the brightenings and velocity episodes are 
visible in our dataset. Figure [4] summarises the radiance 
and velocity evolution in these footpoints, labelled A and 
B in Fig. [3] (TRACE image of 14:42:00 UTC). A and B 
are optimised to obtain a high S/N ratio in the smallest 
region still showing relevant coronal velocities. The area A 
consists of three CDS pixels (6" x 10"), while the area B, 
due to the weak emission of the Fe Xix line in the westward 
ribbon, consists of eleven CDS pixels (224 arcsec 2 ). 



4. Impulsive phase: chromospheric evaporation 

For both footpoints, the radiances of O v and TRACE 
17.1 nm undergo a sharp increase and a rapid decay during 
the impulsive phase. However, Fe xvi and Fe xix present 
a very different behaviour. In A these hotter signatures 
increase and decrease rapidly, while in B Fe xvi becomes 
barely brighter than average and Fe xix begins a slow rise. 
We believe that these differences reflect the different tem- 
perature reached by the evaporated plasma, being hotter 
in B (>10 7 K, cf. §0. 

4.1. Motions in the flaring ribbons 

Footpoint A: Around 14:40 UTC, CDS was scanning 
the east flaring kernel (A), revealing a strong emission 
in O v. Both Fe xvi and Fe xix clearly increased their 
emission with respect to the background. In the 3 pixels 
defining the kernel A downflows of about 20 km s _1 are 
measured in O v, but adjacent pixels show both strong 
upflows (—35 km s _1 , westward of A) and downflows 
(60 km s _1 , eastward of A, see Fig. [5]). The O v down- 
flows appear clearly associated with the bright kernel visi- 
ble in TRACE and in t he HASTA Ha image (see Fig. 1 in 
Maltagliati et al. l2006tl . The same velocity pattern is also 
shown by the He I spectra. 

At this time the Fe xix spectrum obtained integrat- 
ing over A shows a blue-shifted emission with a veloc- 
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ity between —170 km s _1 (single component fitting) and 
—250 km s _1 (double component fitting, see Fig. [51 top- 
right). However, given the noise level in the spectrum, it 
is difficult to establish whether a two component fitting is 
justified. We thus can assume as a lower limit the value 
obtained with a single component fit. The same region 
also shows upward motions around —70 km s _1 in the 
Fe XVI line. Due to the broad NIS-1 profiles, we have not 
attempted to separate the strong Fe xvi background com- 
ponent, resulting in a likely underestimate of the upward 
velocity. 

These upward motions of hot plasma seem to originate 
in the region between the flow patterns of opposite sign 
observed in TR, but from an area where O v motions are 
still downward directed (« 20 km s _1 ). Furthermore, if we 
consider a semi-circular loop with radial legs, the projec- 
tion effect would amount to about —2" in the x-direction 
and thus suggest a link between the upward motions in 
Fe xix and the even stronger downflows seen in the TR 
(the target area is on the west hemisphere, the loop direc- 
tion is east-west, and we assume the Fe Xix emission from 
a region «5 000 km above that of O v and He I ). 

Around 14:44 UTC, CDS scanned again this flaring 
kernel (Fig. [6]). The resulting Fe xix profile (bottom panel 
of the figure) shows a strong blue-shifted component dom- 
inating the whole profile. The corresponding velocity is of 
about —200 km s _1 . Upflows of —70 km s _1 are still mea- 
sured in the coronal Fe xvi. The same velocity pattern 
seen in O v and He I during the previous passage is still 
discernible, with the downflows now weaker (i=s 20 km s ) 
and unchanged upflows. 

Finally, strong coronal upflows are still visible in the 
kernel A when CDS scans the area a third time during 
the impulsive phase, around 14:47:30 UTC (-170 km s _1 
with a double component fitting and —100 km s _1 with 
a single component fitting). O v and He I do not show 
any downflow anymore, while the nearby upflow is still 
present. 

Footpoint B: This footpoint is within the DST FOV. 
Starting from around 14:41 UTC, the difference of Ha 
images obtained at both AA=± 60 pm and ± 150 pm 
shows a red-wing emission excess in correspondence of the 
flaring kernels, resulting in a downflow velocity of 4 to 
10 km s- 1 (cf. §EH). 

Around 14:42 UTC this area was rastered over by CDS. 
Velocities in O v and He I indicate only downwards mo- 
tions (w20 km s _1 ) outlining the tiny bright features vis- 
ible at the higher spatial resolution afforded by Ha data. 
Chromosphcric downflows are measured at the same time 
in the same area, and are indicated with thin black con- 
tours in the upper panels of Fig. [7] Upflows of more than 
200 km s" 1 were measured in Fe Xix by averaging over 
the whole footpoint area. The corresponding profile of the 



Fe XIX line is shown in the bottom panel of the same figure. 
The fitting procedure reveals a blue-shifted component of 
the same strength of the stationary component, with a 
speed of —240 km s . A single component fit provides a 
velocity of about —125 km s -1 . 

During the next raster at 14:45:23 UTC, CDS still 
measures strong upward speeds in the Fe xix line 
(Fig. [4j, cospatial and cotemporal with downflows in the 
low chromosphere (Ha). The radiance of the blue-shifted 
component is comparable to that of the stationary com- 
ponent also at this time. At later times, neither the Fe XIX 
nor the Ha data show any relevant motions. Hence in this 
footpoint, the oppositely directed motions measured in 
low-chromosphere and in hot flare lines, are verified to be 
co-spatial and co-temporal and to persist for at least 200 s. 

Contrary to what could be expected, the Fe xvi line 
does not show any upward velocity within the area B. 
However, we note that at the times when strong motions 
are measured in Fe XIX, the radiance of the Fe XVI is prac- 
tically unchanged with respect to its pre-flare values. This 
also suggests that the evaporated plasma is too hot to gen- 
erate a shifted component strong enough to be measurable 
with respect to the Fe xvi background emission. 

Comparing the two footpoints, we notice that the 
Fexix blue-shifted component is clearly dominating over 
the stationary one only in the eastward ribbon (A), pos- 
sibly due to the smaller area used to obtain reliable line 
profiles. This fact may then support the assumption that 
multiple thin loops are present in the flaring ribbons at any 
time, and that they are heated succe ssively, as in the sim- 
ulations of Warren fe Doschek ( 2005f ). This is further con- 
firmed by TRACE and Ha images that show very clearly 
the existence of multiple footpoints within the flaring rib- 
bons, brightening at different times. Note also that during 
the gradual phase (see § [5]) post-flare loops of 1" to 2" 
diameter clearly emerge as independent features in the 
TRACE images. 

Finally, for both footpoints, the large upflows mea- 
sured in the Fe xix line are spatially related to patches 
of downflows in the TR, themselves closely associated to 
the flaring kernels. However, a significant upflow pattern 
lasting at least 600 s is also present nearby footpoint A, 
but we are at lack of an obvious explanation for it. 



4.2. Momentum balance 

Hydrodynamical simulations of explosive chromospheric 
evaporation predict the equality of momenta between the 
hot plasma moving upward and the dense cold plasma 
moving downward d uring the impulsive phase of a flare 
([Fisher et al.lll985bl ) . If the oppositely directed flows that 
we measure during the impulsive phase indeed signal 
chromospheric evaporation, we should be able to verify 
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Fig. 5. Eastward ribbon (A) at 
14:40:45 UTC. The top left panel 
shows a TRACE image in the 17.1 nm 
band. The bottom left panel shows a 
map of the O v LOS velocity within 
the area encompassed by the black 
box on the TRACE image. On both 
panels isocontours of the O v radiances 
at a level of 5 W m 2 sr _1 are shown 
with white thick lines while the thick 
black isocontours on the velocity 
map indicate Fe xix radiances at a 
level of 90 mW m 2 sr" 1 . The Fe xix 
profile shown on the top-right panel 
is obtained by integrating over three 
pixels along the CDS slit, centred at 
the location indicated by a white x on 
the velocity map (area A defined in 
Fig. [3]). The single pixel O v profiles 
on the bottom right panel refer to 
the locations indicated by a black + 
(thick histogram) and a white + (thin 
histogram) on the velocity map. 
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Fig. 6. Eastward ribbon (A) at 
14:44 UTC. The top left panel shows 
a TRACE image in the 17.1 nm band 
while the Ha image closest in time 
is shown on the top right panel. On 
both images the thin black contours 
indicate chromospheric downflows of 
4 km s _1 . The Fe xix line profile shown 
on the bottom panel is obtained by 
integrating over three pixels along the 
CDS slit (area A defined in Fig. EJ). 



such equality. Earlier attempts h ave been per f ormed , us- ICanfield et al.l (|1990a[) . These authors find agreement be- 
ing BCS/SMM coronal data, by IZarro et al. ( 19881 ) and tween momenta of the oppositely directed plasma (inte- 
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grated over the whole impulsive phase) within one order 
of magnitude, although the lack of spatial resolution in 
coronal data leaves open the issue of co-spatiality of the 
flows. We have shown that chromospheric downflows and 
coronal upflows are co-spatial and co-temporal at least in 
the case of footpoint B, so we can compare the "instanta- 
neous" momenta at 14:42 UTC in this footpoint. 

The calculation of the downward chromospheric 
(Pdown) and the upward coronal (Pup) momenta is de- 
scribed in Appendix [XJ where the uncertainties of the 
measured or assumed physical parameters are taken 
into account. We also note that the calculations of 
both momenta involve totally independent sets of phys- 
ical quantities. For the downward momentum we obtain 
5.1 x lO^gcms- 1 < P down < 4.6 x 10 20 g cm s" 1 . 
The range of values depends from the pre-flare hydro- 
gen density, the chromospheric velocity and the con- 
densation thickness that we can estimate only within a 
factor of two. For the upward momentum, we obtain 
5.9 x 10 19 g cm s" 1 < P up < 4.9 x 10 20 g cm s"\ 
where the lower and upper limits refer to the velocity and 
density values obtained assuming a single or a double- 
component fit of the Fe Xix line, an iron abundance vary- 
ing between photospheric and coronal values and the min- 
imum and maximum value of the coronal filling factor. 
Possible departures from ionisation equilibrium can also 
contribute to the uncertainty of the result. In principle 
we could also attempt to check the momentum balance 
at 14:45:23 UTC, when we still measure co-spatial coro- 
nal and chromospheric flows on the area analysed above. 
However, at this time the chromospheric density has al- 
ready increased to flare values, difficult to estimate with- 
out an ad hoc flare model, which we cannot obtain for lack 
of suitable spectral signatures. 

The equality within one order of magnitude of the "in- 
stantaneous" chromospheric and coronal momenta sup- 
ports the explosive chromospheric evaporation model, 
even though the flare analyzed is quite small and thus 
unlikely to provide the large energy flux requested by the 
simulations. More in general, we note that even with si- 
multaneous, spatially resolved observations we cannot ob- 
tain a verification of the momentum balance better than 
one order of magnitude due to inherent uncertainties in 
the physical parameters. 



5. Gradual phase 

The gradual phase of the flare starts around 14:48 UTC. 
As clearly visible in Fig. SI during this phase no significant 
upflows in coronal lines are measured (nor corresponding 
chromospheric downflows). This indicates the end of chro- 
mospheric evaporation within the flare. 

RHESSI data were available from 14:50 UTC. Images 
were reconstructed in the 3 to 12 keV band, at 14:50, 



14:52 and 15:13 UTC, with an integration time of 60 s us- 
i ng the Pixon algorithm of the standard RHESSI software 
(|Maltagliati et al.ll2006l ). In Fig. M we show the TRACE 
and Ha images at 14:48 UTC, with overplotted, respec- 
tively, the contours of Fe XIX radiance and of the RHESSI 
image in the 6 to 12 keV band at 14:50 UTC (thick solid 
white line). The RHESSI emission encompasses the west- 
ern ribbon and extends eastward outlining a loop struc- 
ture that is also partially visible from the Fe Xix radi- 
ance isocontours (note that the loop top lies outside the 
CDS FOV), but no emission is registered on the eastern 
footpoint (A). Thus, also RHESSI data outline a different 
behaviour of the two flaring ribbons. 

The corresponding integrated RHESSI spectrum 
mainly describes the emission of an area closely overlap- 
ping kernel B. This spectrum (Fig. [51 left panel) reveals 
the broadened emission line feature around 6.6 keV, due 
to lines of highly ionised iron (Fe XXIV - xxvi) observ- 
able in R HESSI data at temperatures above ss 10 7 K 
(jPhillipd l2004h . Indeed, the best fit to the spectrum in- 
dicates a thermal component of T w 10 7 K, as well as a 
non-thermal one with a spectral index 7 sa 10. 7 values 
of the same order have already been obser ved in RHESSI 



data during the decay phas es of small flares (jKrucker et al 



2002; Ha nnah et al.l I2004T) . The presence of non-thermal 



electrons at these low energies, however, is not obviously 
linked to any evaporation process, as no upflows are mea- 
sured in the coronal signatures (bottom panel of Fig. [8]) . 
Around 23 min. later (Fig. [9j right panel) the spectrum 
exhibits only a thermal component with a temperature 
w 8 x 10 6 K. 

A coronal plasma heated above 10 7 K in kernel B might 
explain the very different behaviour displayed by the two 
footpoints, both in radiance and velocity. In A, after the 
end of the impulsive phase, both coronal signatures un- 
dergo a rapid decay, while in B the radiance of Fe xix 
continues its increase up to 15:05 UTC, i.e., well into the 
gradual phase. The Fe xvi shows a similar trend, but with 
a more significant delay, reaching its maximum radiance in 
B almost one hour after the GOES peak. This seems con- 
sistent with the rapid cooling of "hot" plasma (> 10 7 K), 
evaporated during the impulsive phase, to Fe XIX temper- 
atures, and successively, to values characteristic of Fe xvi 
emission. 

In this same footpoint, downward directed motions are 
measured in both Fe xix and Fe XVI around the time of 
their peak emission, although the amplitudes are close to 
the error level (for Fe xix, (20 to 70 ± 50) km s -1 for a 
double-component fitting and (6 ± 12) km s _1 for a single- 
component fit; for Fe xvi, (20 ± 12) km s _1 ). Downflows 
of « (20± 9) km s _1 , are further observed in O v to outline 
the B footpoint around the time of maximum O v emis- 
sion during the gradual phase («15:42 UTC). No obvious 
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Fig. 7. Westward ribbon (B) at 
14:42 UTC. Panels are as in Fig. |1 On 
both images the thin black contours 
indicate chromospheric downflows of 
4 km s _1 and 7 km s _1 . The Fe xix 
line profile shown on the bottom panel 
is obtained by integrating over the 
area within white solid lines in the Ha 
image (area B defined in Fig. [3}. A 
single component fit (plus background) 
would result in an upflow speed of 
(-125 ± 50) km s" 1 
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Fig. 8. Gradual phase at a time when 
the CDS slit was crossing the western 
footpoint. The top left panel show a 
TRACE image in the 17.1 nm band 
while the Ha image closest in time 
is shown on the top right panel. On 
the TRACE image isocontours of the 
Fe xix radiance at levels of (0.11, 
0.22) W m~ 2 sr" 1 are indicated by 
thick white solid lines. On the Ha im- 
age, instead, the white thick contours 
refer to a RHESSI image obtained be- 
tween 14:50 and 14:51 UTC in the 
6 to 12 keV band. The contours are 
traced at levels of 40 % and 75 % of 
the maximum. The Fe xix line profile 
shown on the bottom panel is obtained 
by integrating over the areas within 
white solid thin lines in the Ha image 
(area B). In this case the radiance of 
the rest component (dashed line) was 
set to zero. The dotted line indicate the 
background (lines plus constant back- 
ground) defined from pre-flare spectra. 
The thick solid line indicate the fitted 
profile. 



downflows are instead measured in A during the gradual Finally, combining gradual phase RHESSI data with 
phase. CDS measurements up to 14:50 UTC, we evaluate the 



Teriaca et al.: Dynamics and evolution of an eruptive flare 



11 



contribution of the evaporated chromospheric plasma in 
filling the flare loops. Using the Fe xix data, we can esti- 
mate the total number of electrons transported with evap- 
oration during the impulsive phase, assuming a symmetric 
evaporation on the two footpoints: 



N U p — n e 2 S c , 



f v cor At, 



(1) 



where n c is the electron density of the upflowing plasma, 
/ the coronal filling factor, 2 S cor the area of the ob- 
served upflows (both footpoints), v cor the velocity, and 
At the duration of the evaporation. The electron den- 
sity and the upward velocity are the same as in § 14.2 
(see also Eq. IA.3|) and depend on whether we use one 
(v C or — —125 km s _1 ) or two (v cor = —240 km s _1 ) com- 
ponents to fit the Fe xix line profile at 14:42:03 UTC and 
upon the assumed iron abundance. We assume such densi- 
ties and velocities as typical throughout the evaporation 



2S C 



2 x 1.2 x 10 18 cm 2 accounts for the fact that the 



plasma is flowing from both footpoints. Finally, At — 300 s 
is taken equal to the duration of chromospheric downflows, 
for which we have a better temporal resolution. 

With these values, for the two cases of single and 
double fitting and for both iron abundances of 7.5 and 
8.1, we obtain the minimum and maximum estimate of 
N up = [3.8,10.] xl0 37 \/7 electrons. 

The above values can be compared with the total num- 
ber of electrons in the entire loop system after the end of 
the evaporation phase. This can be written as: 



N, 



'EM 



loop 



RHESSI 



f S cor L 



loop 



f S C or L 



loop J 



(2) 



0.8 x 10 8 cm 6 is the emis- 



where EM RHE ssi 
sion measure obtained from the RHESSI spectrum at 
15:13 UTC (only thermal component, see Fig. [9]). From 
the TRACE images around 16:20 UTC a loop length 
Lioop of « 9.2 x 10 9 cm can be derived, from which 
Nioop = 9.4 x 10 37 \/7 electrons. Since both N up and Ni oop 
have the same dependence on y/J, the evaporated material 
can account for 40 % up to 100 % of the material inside 
the loop system. This implies that the evaporated plasma 
largely contribute to the density of the post-flare loops. 



6. Cooling of post-flare loops 

The whole evolution of the flare outlines the filling of mag- 
netic loops by plasma evaporated from the chromosphere 
during the impulsive phase. Such loops become visible in 
all of our signatures during the extended gradual phase, 
with progressively cooler features appearing at later times 
(see Fig. Fig.[lI3reports the typical temperatures corre- 
sponding to our signatures, assumed as those correspond- 
ing to the maximum of the contribution function, vs. the 



peak time of the loop radiances. Moreover, the tempera- 
ture of the thermal spectrum used to fit RHESSI data is 
also reported at the two times of acquisition. 

We compare here the temporal evolution defined by 
the data with the results of the analytical formulae pro- 
vided by I Car gill et al. () 19951 ) for the cooling time of a flare 
plasma. The first RHESSI spectrum at 14:50 UTC still 
shows a non-thermal component (see § [5]) , suggesting fur- 
ther energy input at the beginning of the gradual phase. 
Hence we begin the analysis of the cooling at 15:13 UTC, 
when RHESSI shows only a thermal component. The first 
data point corresponds to a T » 8 x 10 6 K, and an av- 
erage electron density of 8.5 x 10 9 cm -3 (for / = 1). We 
used the formulae provided for the static case since after 
15:13 UTC no evaporative motions are visible anymore. 
The calculations were performed following two different 
approaches. 

In the first case a constant density was assumed 
throughout the entire co oling period and th e complete 
radiative loss function of Rosner et al.1 ( 19781 ) was used 
in the calculations . For this purpose Eqs. 7S and 8S of 
Cargill et al. (1995) have been generalised to any value of 



the exponent a of the power law. Using such equations, 
we computed the time until when conduction dominates 
over radiation. The temperature evolution of the cool- 



ing plasma was obtained through Eq. 3a of I Car gill et al 



([19951 ) before this time and Eq. 4 after, with the appro- 
pr iate values fo r each piece of the radiative loss function 
of iRosner etal] (|l978f) . The thin lines in Fig. [TDI show the 
results of these calculations for three values of the filling 
factor / (0.1, dot-dashed; 0.3, dashed; 0.5, dotted). The 
solid part of the curves indicates the times when conduc- 
tion dominates, while the non-solid part indicates the time 
when the plasma cools radiatively. 

In the second approach we assumed that the den- 
sity is constant only during the c onductive phase while 



n oc \/T in the radiative phase ( Serio et all Il99ll ). We 
also assumed the simplified radiative loss function P ra d = 
1.7 x lQ- 19 n 2 T- x / 2 that fits, within a factor two, the 
CHIANTI radiative loss function for coronal abundances 
for 5 < log(T/K) < 7. During the time dominated by 
conduction, the calculations are similar to the previous 
ca se. After this time, being a — —0.5, we can use Eq. 6 
of ICargill et al.l (|1995l ) to calculate the temperature evo- 
lution during the time when radiative cooling dominates. 
The results are represented by the thick curves in Fig. [TQl 
for the same values of /. Again, the solid part represent 
the conductive cooling while the non-solid part is the ra- 
diative one. 

Both approaches show that the observed curve can 
be explained by the cooling of the hot flare plasma in- 
jected into the loop system only if a volumetric fill- 
ing factor of w 0.2 to 0.5 is assumed, a value well in 
agreement with the estimated filling factors of cooling 
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Fig. 9. RHESSI spectra in the 4 to 10 keV energy range at 14:50 UTC (left panel) and at 15:13 UTC (right panel). 
Notice the presence of a broadened emission feature at about 6.6 keV corresponding to a group of emission lines due 
to Fe xxiv - xxvi. 



loops given by lAschwanden et al" (2003). Wc notice that 
the estimated cooling time in the temperature range 
4 < log(T/K) < 5 is much shorter than the observed 
one. This is expected because our initial assumptions are 
no longer valid in this range, in particular due to recom- 
bination of hydrogen. 

7. Conclusions 

We presented comprehensive observations of a small erup- 
tive flare (GOES class C2.3) that developed on 2002 
August 11 in NOAA region 10061. A wide array of instru- 
ments and diagnostics allowed us to follow the evolution 
of the flare from the chromosphere to the corona during 
both the impulsive phase and the long gradual phase. The 
eruption of a portion of a filament, located on the appar- 
ent magnetic inversion line, was signalled by an upward 
velocity of about —30 km s _1 detected at 14:36:33 UTC 
by CDS in the O V line. A few minutes later the foot- 
points of a large loop system, subsequently visible also in 
the TRACE images, began to brighten in all signatures, 
reaching their maximum brightness around 14:42 UTC. 

During the impulsive phase of the flare, large coronal 
upflows of —200 km s _1 or more are measured from the 
Doppler shift of the Fe xix line profiles on the ribbons at 
the footpoints of the flaring loop system. Interestingly, the 
two footpoints show quite a different behaviour. Footpoint 
A shows a fast (within our time resolution) rise of the 
radiances of all our signatures, followed by a slower decay 
of the Fe xvi and Fe Xix emission. On this footpoint, 
during the impulsive phase, we can select a small area 



(3 CDS pixels) where a blue-shifted component clearly 
dominates the Fe XIX profile. Upward motions of about 
— 70 km s _1 (likely underestimated) are also measured 
from the Fe xvi CDS data. These coronal upward motions 
occur in the location where downflows are observed in the 
O v and He I lines for at least 200 s but, unfortunately, 
no Ha Dopplergrams are available for this footpoint. 



On footpoint B instead, the radiances of Fe xix and 
Fe xvi lines increase very slowly and reach their maxima 
respectively 20 min and 60 min after the start of the flare. 
RHESSI data at 14:50 UTC show the emission in the 3 
to 12 keV band to be mainly associated to this footpoint, 
with the spectrum indicating a temperature of the thermal 
component of 10.4 MK, i.e. larger than the Fe xix forma- 
tion temperature (8 MK). This would explain the slow rise 
of Fe XIX and Fe xvi radiances as due to progressive cool- 
ing of this "hot" plasma. During the impulsive phase, the 
blue-shifted and the stationary component of the Fe xix 
line profile on this area are of comparable strength, and 
no clear motion is measured from the Fe xvi line. The O v 
and He I lines show downward (« 20 km s _1 ) motions out- 
lining the flaring kernels. For this footpoint, within the 
FOV of the UBF imager, we observed simultaneous Ha 
downflows that we estimated to be about (7± 3) km s -1 . 
To our knowledge, this is the first time that co-spatial, co- 
temporal, oppositely directed flows during the impulsive 
phase of a flare are measured in low chromosphere and 
corona, supporting the model of explosive chromospheric 
evaporation. 
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Fig. 10. Cooling curve obtained from the time of maxi- 
mum appearance of progressively colder features. RHESSI 
temperature at 14:50 UTC refers to the thermal compo- 
nent (see left panel of Fig. [5]). Each set of curves (thin 
and thick) refer to three different values of the volumetric 
filling factor / (0.1, dot-dashed; 0.3, dashed; 0.5, dotted). 
The solid part of the curves indicates the times when con- 
duction dominates, while the non-solid part indicates the 
time when the plasma cools radiatively. The thick curves 
are obtained assuming constant density only during the 
conductive phase, while the thin curves are obtained as- 
suming constant density through the entire cooling (see 
text). 



Such direct measurement gave us the possibility to 
compare the "instantaneous" momenta of the coronal and 
chromospheric moving plasma, whose equality is required 
in the hypothesis of explosive evaporation. We found that 
such momenta a re equal within one order of magnitude . 
Previous works ( Zarro et al. 19881 : Canfield et al. 1990a ) 
verified the equality of the momenta within the same lim- 
its using data without spatial resolution and integrating 
over the whole impulsive phase. Although we can utilise 
spatially resolved data, it remains very difficult to obtain 
more accurate estimates of the momentum balance, due to 
the inherent uncertainties in physical parameters such as 
coronal electron density and abundances, chromospheric 
hydrogen density and condensation height. To improve the 
knowledge of these parameters a larger number of spectral 
lines both in the EUV (including density diagnostics) and 
in visible/near-infrared (to build a semi-empirical model 
of the chromosphere) are necessary. However, this is likely 
to jeopardise the temporal and spatial resolution needed 
to study the impulsive phase of a flare. New space instru- 
ments with larger effective areas, such as the forthcoming 



EIS on Solar B, should improve the situation on the EUV 
side. 

A further element in support of the chromospheric 
evaporation was given by RHESSI data acquired during 
the gradual phase, from which we obtained an estimate 
of the mass in the coronal loop system. Combining the 
velocity and radiance data from CDS, we could also esti- 
mate the mass evaporated during the impulsive phase of 
the flare, and found that it provided a large contribution 
to the density of the X-ray emitting coronal material. 

Finally, during the extended gradual phase, the mag- 
netic loops filled by evaporated plasma become visible in 
all of our spectral lines and band-passes, with progres- 
sively cooler features appearing at later times. The ob- 
served temporal evolu tion can be r e produ ced using the 
analytical formulae of ICargill et al. (1995). Initially the 
cooling is dominated by conduction, while radiative losses 
take over later on. A cooling time comparable to the ob- 
served one can be obtained only assuming a volumetric 
filling factor of as 0.2 to 0.5. This in turn is consistent 
with the the presence of multiple thin loops heated suc- 
cessively, suggested by the observed dominant blue-shifted 
component in the Fe XIX line profile only when we can use 
a small area to obtain reliable data. 

To end this paper, we would like to make some remarks 
about the complexity of the TR behaviour. While our ob- 
servations show downflows in the flaring kernels also from 
TR lines, in qualitative agreement with so me hydrody- 



namical simulations of explosive evaporation ([Fisher et al 



1985al fbt). we note that rece nt hydrodynamic m odels of 



chromospheric evaporation ( Allred et al. 2005) always 
predict upflows for the TR plasma. Several observations 
obtained with CDS during the impulsive phase of flares, 
with programs privilegeing either the spatial or the tempo- 
ral resolution, offe r a wide range of res ults. Upflows have 



been reported by iTeriaca et al.1 (|2003l ). using an obser- 



vational sequenc e simil ar to the one presented in this pa- 
per, and lBrosiud (|2003l ). that adopted a very high cadence 
program, but with a narrow field of view and limited spa- 
tial resolution. Upflow is also measured in our data dur- 
ing the impulsive phase, in the eastern footpoint nearby 
the area A but not appear to be obviously connected to 
the coronal upflows.. Do wnflows are instead reported by 
iBrosius fe Philli pj ( 20041 ) . using the same obse rving pro- 
gram as lBrosiud (|2003L and lKamio et al.l(|2005h that anal- 
yse four flares in a larger FOV but with intermediate tem- 
poral resolution. Although at much lower spatial and tem- 
poral resolution, the same contradictory indications were 
obtained wi th SMM and OSO-8 instruments (see the dis- 



cussion m 



Teriaca et al.l 120031 ) . It is likely that the be- 
haviour of the TR is strongly related to the structure of 
the atmosphere depending upon the local magnetic topol- 
ogy. Space observations with substantially higher spatial 
resolution than now (< 1") are most likely necessary to 
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better understanding the dynamics of the TR and provide 
constraints to theoretical models. Further, we would like 
to point out that hydrodynamical simulations are gener- 
ally computed for energetic flares (GOES class X-M) with 
a much harder spectrum for the impinging electrons (7 
from 4 to 6) than normally observed with RHESSI (7 from 
7 to 10) in small flares such as the one studied here. 

Acknowledgements. The authors thank Prof. E. Priest for sug- 
gesting calculating the amount of plasma evaporated into the 
loop and Dr. K. Wilhelm for carefully reading the manuscript. 
We are also grateful to Dr. V. Andretta and Dr. U. Schiihle for 
useful comments and suggestions. Many thanks are due to the 
NSO and CDS staffs for their help and assistance in acquir- 
ing the data here analysed. We would also like to thank the 
RHESSI staff for their help in analysing the data. SOHO is a 
mission of international cooperation between ESA and NASA. 
NSO is operated by AURA, Inc., under cooperative agreement 
with the NSF. 

Appendix A: Momenta calculation. 

We give here a detailed description of equations and phys- 
ical parameters used to calculate the downward and up- 



with n c given by 



ward momenta P<i 



and P u 



The mass density can be written as p = fim p tit, where 
11 is the mean particle weight, m p = 1.67 x 10~ 24 g is 
the proton mass, and tit the total particle number den- 
sity. For a fully-ionised plasma of hydrogen with 10 % 
helium, 11 — 0.61 and tit = E91 n e , while for a neutral 
gas, ti = 1.27 and tit = IT nu, where nu is the hydro- 
gen number density. Downflow momentum in the chromo- 
sphere (neutral gas) can then be written as: 

Pdown = 1-27 m p 1.1 n ch S ch v ch Ah f ch , (AT) 

where n c h is the pre-flare chromospheric hydrogen density, 
ass umed to be (2.2 to 4 0) x 10 13 cm -3 from model VAL 



F ( Vernazza et al 1 ll98ll) . The with detectable 

chromospheric downflows is 2.5 x 10 17 cm 2 . Within this 
area, the average downflow velocity, v c h, is about (7 ± 



3) km 



A/i represents the thickness of the chromo- 



spheric condensation. We assume that A h between a 
lower limit of sa 100 km , given by dynamic simulations 
( Abbett fe Hawlevlll999l) and an upper limit of 200 km, 
obtained in semi-empirical models of a small flar e that in- 
cludes the velocity fields ( Falchi fe Mauad 2002). Finally, 
f c h is the chromospheric filling factor. Since the spa- 
tial scale of the Hcv images is 0.5" per pixel, smaller 
than the size of the velocity patches in the chromosphere 
(2" to 3"), we assume f c h = 1. With these values, we ob- 
tain 5.1 x lO^gcms- 1 < P down < 4.6 x 10 20 g cm s" 1 . 
The lower and upper limits refer to the minimum and max- 
imum values for the hydrogen density (n c h), v c h, and A h. 

At coronal level, the momentum of the upflowing 
plasma is given by: 



EM 



f S c , 



(A.3) 



In the above equations S cor is the area of the re- 
gion from where blue-shifted profiles arise, L the thick- 
ness of the emitting region, / the coronal filling fac- 
tor, v cor the plasma speed at coronal level, and EM 
the emission measure computed from the observed ra- 
diance using the CHIANTI database ijDere et all Il997 : 
Landi et al.lll 999) wit h ioni sation equilibrium calculations 



from lMazzotta et all {l998) 



The radiance of the Fe xix 59.2 nm line is proportional 
to n e 2 for 8 < log n e /cm~ 3 < 12 and it is, hence, suit- 
able for emission measure calculation. The iron abundance 
Arpe] is 8.1 in the corona (jFeldman et al.lll992h . However, 



because the plasma is ablated fr om the chromosphere, the 
photospheric abundance of 7.6 ( Grevesse fe Sauval 19981 ) 
may also apply. Thus, we will assume an abundance vary- 
ing between 7.5 and 8.1. 

Coronal upflows are measured in eleven effective CDS 
pixels, so S cor — 1.2 x 10 18 cm 2 . L is taken to vary be- 
tween 1.1 xlO 9 cm (= \/S cor ) and 2xl0 9 cm (obtained 
estimating the geometrical depth of the flaring loop along 
the line of sight, assuming the loop is perpendicular to the 
solar surface). The electron density, n c , of the evaporat- 
ing plasma (from the emission measure) and the coronal 
velocity, v cor , depend on whether the Fe XIX spectra are 
fitted with one or two components. The computed values 
of n c are between 3 x 10 9 and 8 x 10 9 cm -3 for f=l. 
Finally, we note that in the extreme case of magnetic flux 
tubes not expanding with height, the ratio S c h/S cor yields 
a lower limit of / = 0.2, so we assume 0.2 < / < 1. 

For both values of /, P up is calculated with the 
values of EM and v cor obtained in both cases of sin- 
gle and double component fitting and for both es- 



timates of L and A 



[Fe]- 



The 



maximum 



and 



Pup = 0.61 m p 1.91 n c S c 



L f v c 



(A.2) 



mum of the sixteen obtained values are taken to give 
5.9 x 10 19 g cm s" 1 < P up < 4.9 x 10 20 g cm s" 1 . 
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